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Abstract 

The intrinsic spin-Hall effects (SHE) in p-doped semiconductors [S. Murakami et al, Science 
301, 1348 (2003)] and two-dimensional electron gases with Rashba spin-orbit coupling [J. Sinova 
et al, Phys. Rev. Lett. 92, 126603 (2004)] have been the subject of many theoretical studies, 
but their driving mechanisms have yet to be described in a unified manner. The former effect 
arises from the adiabatic topological curvature of momentum space, from which holes acquire a 
spin-dependent anomalous velocity. The SHE in the Rashba system, on the other hand, results 
from the momentum-dependent spin dynamics in the presence of an external electric field. The 
two effects clearly appear to originate from distinct mechanisms. Our motivation for this article is 
to address this apparent disparity and, in particular, to seek a unifying description of the effects. 
In this endeavor, we consider the explicit time-dependence of SHE systems starting with a general 
spin-orbit model. We find that by performing a gauge transformation of the general model with 
respect to time, a well-defined gauge field appears in time space which has the physical significance 
of an effective magnetic field. This magnetic field is shown to precisely account for the SHE in the 
Rashba system in the adiabatic limit. Remarkably, by carefully analyzing the equation of motion 
of the general model, this field component is also found to be the origin of the anomalous velocity 
due to the momentum space curvature. Our study therefore unifies the two seemingly disparate 
intrinsic SHEs under a common adiabatic framework. 

PACS numbers: 03.65.Vf, 03.65.-w,73.63.-b 
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I. INTRODUCTION 



The spin-Hall effects (SHE) are a family of phenomena in condensed matter systems 
in which an applied longitudinal electric field gives rise to a transverse spin current. The 
spin current arises from the transverse separation of spin species in the system, and the 
physics driving the separation mechanism can be rather distinct across different systems. 
The earliest prediction of transverse spin separation was made by D'yakonov and Perel' in 
the 1970's,- who studied the spin-dependent scattering mechanisms of carriers with localized 
impurities. The class of SHEs which occur as a result of the spin-orbit interaction of carriers 
with impurities are referred to in the literature as extrinsic. Conversely, there are intrinsic 
forms of the SHE, which have become an active field of research in more recent years,- 1 ^ 1 ^ 
following two seminal papers: the first by Murakami et al. which predicts a SHE of holes 
in p-doped semiconductors described in the Luttinger model,- and the other by Sinova 
et al. in two-dimensional electron gases (2DEG) formed in semiconductor heterostructures 
with Rashba spin-orbit coupling. 7 Both effects are finite in the absence of disorder, and are 
characterized by a pure transverse spin current generated from the spin-orbit coupling (SOC) 
in the band structure of the system. Previous studies^ have shown that the SHE described 
in Ref. for infinite Rashba systems vanishes when one includes vertex corrections to model 
the effects of impurity scattering. However, the effect may still be manifested in finite-sized 
systems^^ 1 ^ such as in mesoscopically confined 2DEGs [lOj (e.g. in quasi one- dimensional 
quantum wires), or in the presence of magnetic impurities.— In this work, we will analyze 
the SHE in Rashba systems without considering vertex corrections (i.e. in line with the 

n 

original treatment in Ref. [7[]). Instead, we will provide a phenomenological explanation of 
the vanishing SHE based on our own analysis later in the paper. In constrast, the SHE in 
the Luttinger hole system is robust to vertex corrections.— 



It is intriguing that the two SHEs in Refs. 



6| and although both intrinsic in nature, 



appear to originate from distinct mechanisms. The former effect in Ref. p] is an adiabatic 
effect described via a gauge potential which arises from the relaxation of hole spins to an 
effective magnetic field in momentum (k) space. The topological Berry curvature^ of the 
gauge potential has the physical significance of a magnetic field in momentum space, and 
affects the trajectory of carriers in much the same way as a classical magnetic field does in 
real space. Here, the resulting magnetic Lorentz force in momentum space manifests itself 
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as an additional (anomalous) velocity in real space. The semiclassical equations of motion 
of carriers in the presence of the Berry curvature have been derived previously— and will be 
revisited for the Luttinger system in this article. It is found that the real space trajectory 
of holes along the transverse direction is spin-dependent, thus resulting in a finite SHE. The 
SHE in Rashba systems^ 7 on the other hand, was derived originally from a semiclassical 
analysis of electron spin dynamics in a Rashba 2DEG system, with no apparent relation to 
the fc-space topology. In the analysis,- electrons were found to gain a momentum-dependent 
out-of-plane spin polarization in the presence of Rashba SOC and an external electric field, 
leading to a transverse separation of spins. It is often stated in the literature that the effect 
arises from the ^-anisotropic precession of spins.- 1 ^ 1 ^ 1 ^ As part of the motivation for this 
paper, we will clarify the mechanism and show that the effect is in fact an adiabatic effect 
in which spins become aligned to momentum-dependent effective magnetic fields. 

From a heuristic viewpoint, the physical mechanism of the two SHEs ^J3] are clearly 
distinct: in the former, carriers acquire a spin-dependent anomalous velocity, whilst in the 
latter they acquire a momentum-dependent spin-polarization. Our motivation for this paper 
is two-fold. Firstly, we ask whether the SHE in Rashba systems can also be formulated within 
an adiabatic framework, and, secondly, whether the physical mechanisms of the two SHEs 
can be unified. In order to describe the Rashba SHE under an adiabatic formulation, it is 
instructive to make note of several points: (1) the adiabatic Berry curvature of momentum 
space of the Rashba system vanishes except as a 5-function singularity at k = 0.— Therefore, 
the spin-dependent anomalous velocity in the Rashba system vanishes for electrons with 
k 7^ and thus does not contribute any transverse spin current. In contrast, the spin-Hall 
current in the Luttinger system results from the Dirac monopole curvature (~ k/\k\ 3 ) of 
momentum space, (2) in the Rashba SHE, spins become tilted out-of-plane which appears 
contradictory to the adiabatic regime whereby they are assumed to follow perfectly the in- 
plane Rashba field, and (3) although the Berry curvature in the Rashba system exists only 
at a singular point in fc-space, the resulting Berry phase is finite, and previous studies have 
shown that the spin-Hall conductivity in the Rashba system is related to the Berry phase 
through the Kubo formula.- 1 ^ 

In this article, we find that the above remarks (l)-(3) can be consolidated into a consistent 
adiabatic theory that emerges from a consideration of the explicit time-dependence of SHE 
systems. In particular, a gauge field Ao(t) naturally appears in time space upon applying 
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a unitary transformation to the system, which has the physical significance of a magnetic 
field in the transformed system. This magnetic field couples to the electron spin, and is 
shown to precisely account for the SHE in the Rashba system in the adiabatic limit; here, 
this limit amounts to the spins following the direction of the sum of the Rashba field and 
the new effective field. Furthermore, the Berry phase can be equivalently expressed in terms 
of the adiabatic components of the gauge field Ao(t). Thus, a gauge field description can 
be attributed to both intrinsic SHEs, although their respective gauge fields are defined in 
different spaces (momentum and time). 

Having identified that both intrinsic SHEs arise from gauge fields in the adiabatic limit, 
we finally embark on the problem of unifying the physical origin of the two effects. In 
the presence of an external electric field, the momentum and time spaces become coupled 
through the usual drift equation of charged carriers. Remarkably, by analyzing carefully the 
equations of motion of a general SOC model, it is found that the anomalous velocity due to 
the Berry curvature in momentum space is in fact a direct result of the effective magnetic 
field component arising from Ao(t). In this sense, the common origin of the two seemingly 
disparate SHEs is clarified. 



II. THEORY 

A. Carrier dynamics in the presence of Berry's curvature in momentum space in 
spin-orbit coupling systems 

1. Holes in the Luttinger system 

Let us briefly review the mechanism for the SHE of holes in p-doped semiconductors 
reported in Ref. 6j. The effective Luttinger Hamiltonian for holes in the valence band of 
conventional semiconductors is given by— 

^Lutt. = y (71 + |y) - 1$ ■ Sf + V(r), (1) 

where 71, 7 are valence-band parameters defining the effective hole masses, k is the momen- 
tum operator, S is the vector of spin-3/2 matrices, and V = V(f) is the potential energy 
(in our notation, a hat (6) signifies an operator while an over-arrow (v) signifies a vector). 
The holes described by (OP) have a well-defined chirality, A = h~ l k ■ S/\k\. Because of the 



chirality-squared term in the Hamiltonian, states with opposite signs of chirality (A = ±1/2 
and A = ±3/2) are degenerate (they correspond to the light-hole and heavy- hole bands, 
respectively). In the presence of an external electric field E, the potential energy term is 
V(r) = eE ■ f , where — e is the electron charge. Parameterizing the momentum vector 
k = | A; | (sin ^ cos 0, sin # sin 0, cos6), we proceed to define a 4 x 4 unitary matrix U(k), 

U(k) = exp (i9(k)S y ^ ■ exp (i(()(k)S^ , (2) 

which aligns the reference spin axis to be along the direction of k, i.e. it satisfies the 
relation U(k) (k-SjW(k) = \k\S z . The effective, diagonalized Hamiltonian Ti,' Lutt = 
U{k)'H^ u x,\,.U^ (k) then reads 

nU. = y (71 ± ~t) - l\k?S 2 z + U{k)V(f)U\k). (3) 

In the last term of ([3]), the position operator f = id k acts as a partial derivative in momentum 
space, and we obtain from the fc-dependence of U : 

U(k) (eE ■ id k ^ U\k) = eE ■ (r + iU$)d k U\k)\ . (4) 

Thus, under the local transformation, the position operator transforms into covariant form, 
f — > R = f — A(k), where A(k) = —iU{k)diJJ^(k) is a gauge field in reciprocal space. 
Thus far, the transformation 7iLutt. — > ^Lutt. i s exact. Being a pure gauge field, A(k) 
induced by the transformation has no associated curvature. Assuming adiabatic transport, 
in which we neglect mixing between the light-hole and heavy-hole bands, and applying an 
Abelian approximation within each hole band, we are left with only the diagonal gauge field 
components of the respective 2x2 hole band subspaces. Explicitly, the Abelian gauge fields 
are given by 

A &d -(k) = -Acos#V k 0, (5) 

where the superscript ad. denotes adiabatic transport, and A is the hole chirality. The 
corresponding gauge invariant quantity (which is thus related to a real physical effect) is the 
curvature tensor Q(k), defined by 

n k (k) = d kt At;(k)-d k] At(k). (6) 

— # 

The curvature il(k) above is frequently called the Berry curvature in momentum space. In 
the present simple calculation reveals that the Berry curvature is 

m = x±-, (?) 
\k\ 3 



i.e. it is a Dirac monopole with strength eg = X. It turns out that the /c-space curvature 
(EJ) has important implications on carrier dynamics. In particular, Q(k) can be regarded as 
a magnetic field in fc-space, which gives rise to a fc-space Lorentz-type force. The modified 
semiclassical equations of motion for carriers in the presence of a non-trivial curvature in 
fc-space have been derived elsewhere to be:— 

hk = -eE, (8) 
f = ^V k e- k x Q(k), (9) 

where the over-dot signifies time differentiation, and e is the energy eigenvalue of the system. 
The final term in ([9]) is the Lorentz-type force in fc-space, and is equivalent to an additional 
velocity of electrons corresponding to the so-called anomalous Karplus-Luttinger term.— 
Substituting the expression for the curvature ([7j) into the equation of motion, the anomalous 
velocity component is given by 

t» anom - = -Xkx 4-, (10) 
|Jfe|3 

which is perpendicular to both the applied electric field E and k. Since the chirality of the 
holes has sign A > (< 0) for hole spins (anti-)parallel to the electron momentum, the 
anomalous velocity is also perpendicular to the spin S, and points along opposite directions 
depending on the sign of the chirality. This transverse separation of the spins gives rise to 
the SHE of holes in the Luttinger system. 



2. Conduction electrons in the Rashba system 

We now analyze the Berry curvature in momentum space for the case of the linear Rashba 
SOC,— >2i which is present in two-dimensional electron gases formed in semiconductor het- 
erostructures. We begin with the generalized spin-orbit Hamiltonian, 

H=^-W-B(k)+V(r), (11) 

where m is the effective electron mass, 7 is the SOC strength, a = {a{\ is the vector of 
Pauli spin matrices, B{k) is a momentum- dependent effective magnetic field, and V(f) = 
eE ■ r in the presence of an external electric field E. The above Hamiltonian captures the 
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physics of many other types of SOC, including the linear and cubic Dresselhaus^ and strain- 
induced^ 1 ^ SOC systems. The Luttinger Hamiltonian ([!]) can also be transformed to be of 



this general form when re-cast in terms of the SO (5) Clifford algebra as was done in Ref. [28], 
although this representation is in a 5-dimensional space rather than the usual spin-| space. 
The single particle eigenstates of the Hamiltonian are of the form \ip±) = exp (ik ■ r^%±(fc), 
i.e. a product of the spatial plane wave state and the spinor part which encodes the electron 
spin state. For any k, the spin degeneracy is lifted between the two eigenstates \ip±), which 
have corresponding spin-orbit energy eigenvalues of e± = ±'~f\B(k)\. Let us rotate the 
reference spin axis such that it points along the direction of the spin-orbit field B(k), i.e. we 
diagonalize the Hamiltonian with respect to B(k). By parameterizing the spin-orbit field 
in terms of spherical angles, B = \B\(sin9 cos <p, sin 9 sin 0, cos 6*) = \B\n, where 9 and <fi 
are explicit functions of k, the diagonalization may be achieved through the SU(2) rotation 
matrix U = U{k) given by Eq. ([2]) but with the replacements S y — > cr y /2, S z — > a z /2. 
However, the choice of U for the diagonalization is not unique: for convenience we shall 
adopt another rotation matrix given by^ 

U(k) = m(k) ■ a, (12) 

where m = (sin | cos 0, sin | sin 0, cos |). The effective, diagonalized Hamiltonian is given 
by 

U' = UHU ] = |- - Wz\B{k)\ + UV{f)U ] . (13) 

The <j z Pauli matrix in the diagonalized spin-orbit term represents the two spin states either 
parallel (the ground state) or anti-parallel (the e + state) to the spin-orbit field B(k). In the 
last term of ([TBI , the position operator f = idk is transformed into the covariant form of 
Eq. (j3J), i-e. f — > R = r — A(k), where A(k) is a SU(2) gauge field in reciprocal space. From 
Eqs. (HI) and ( fl2|) . the gauge field components can be represented in terms of the fa- vector 
and the Pauli spin matrices, i.e. 

A ki (k) = {m x d ki m) ■ a = A ki ■ a, (14) 

and in terms of the n-vector by replacing in the above equation with 

where i = x,y,z are real space coordinates. Up to this point, the transformation of the 
Hamiltonian is general. We now impose the adiabatic approximation, in which mixing 



between the two eigenstates of the diagonalized Hamiltonian is neglected. Mathematically, 
this corresponds to retaining only the diagonal terms of A(k), i.e. the a z coefficients in Eq. 
(I14p . from which we obtain an Abelian gauge field known as the Berry connection, ^4 ad '(fc, s). 
A ad '(k, s) has two values, representing the two spin states, s = ±1, of the diagonalized 
Hamiltonian (we denote the ground state as s = +1), and which correspond to the diagonal 
terms of Aki{k). Explicitly, the Abelian gauge field is given by 

^•(£, s ) = -£(l-cos0)V k 0. (16) 

The curvature tensor Q(k) of this connection, defined by Eq. ([6]), is invariant with respect 
to the gauge transformation U(k). From the definition, it is clear that f2(fc) respects the 
same symmetry in s as the connection, i.e. Q(k, s) = —Q(k, —s). 

In principle, one can define the curvature Q in any arbitrary space. For example, in the 
special case of the magnetic field space B, the Berry curvature has the classic form of Dirac's 
monopole,— 

— * 

n(B) = s ^-. (17) 

2\B\ 3 

The above relation is general and applies to any SOC system. One can transform the 
curvature from 5-space to any other space (e.g. fc-space) by using the relation:— 1 ^ 

(dB dB \ 
^x— j, (18) 

where is the Levi-Civita symbol. Generally, the curvature in momentum space Q(k) is 
not the Dirac monopole field (although it still is for the case of the Luttinger Hamiltonian). 
The actual form of Q(k) depends on the fc-dependence of the effective magnetic field. We 
saw in Eq. (Q how this curvature Q(k) gave rise to an anomalous velocity which resulted 
in the SHE in p-doped semiconductors. However, the same reasoning cannot be applied to 
the SHE in the Rashba system, as f2(/c) in this system is vanishing (for k 0) as we outline 
below. 

The Hamiltonian in the presence of the Rashba SOC is given 

2 

2m 

where a is the Rashba spin-orbit coupling parameter expressed in units of eVm. The ef- 
fective magnetic field is given by B R {k) = (k y ,—k x ), and the eigenvectors are \k, ±) = 



H R = ^— + a ( k x a y - k y a x ) , (19) 
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1/v^exp yik ■ rj (^fik~ 1 (k x — ik y ), 1) T , with the corresponding energy eigenvalues of e± = 
±ak where k — \k\ — \Jk"l + is the in-plane wave- vector magnitude. In momentum 
space, the effective magnetic field Br is directed along 9 = tt/2 and <fi = tan _1 (— k x /k y ), and 
from Eq. ( fl6l) the Berry connection is given by A^ik, ±) = ±^{-k y , k x , 0). Evidently the 
curvature ([6]) of this connection is trivial, i.e. = over the entire fc-space, except at 

the singularity point at k = where the /^-component of the curvature has non-vanishing 
value ±7r, i.e. the curvature is of the form Q(k, ±) = (0, 0, ±7r5(/c)).— Thus, conduction 
electrons having a finite momentum in the 2DEG plane do not experience any Lorentz-type 
force in /c-space, as is the case for holes in the Luttinger system. Furthermore, even if this 
force existed, it would only separate the Rashba SOC eigenstates (whose spins lie entirely 
in-plane) in the transverse direction, and thus cannot explain the out-of-plane spin polar- 
ization acquired by the electrons in the SHE. Now the question arises as to whether the 
SHE in Rashba systems can be described within a gauge field framework. In particular, the 
out-of-plane spin polarization seems to suggest the presence of an additional magnetic field 
in the system. It turns out that such a gauge formulation does exist, but one must turn to 
another parameter space, namely the time space. 

B. Time component of the gauge field in spin-orbit coupling systems with an 
electric field 

When considering the temporal evolution of a quantum system, the unitary transforma- 
tion is explicitly time-dependent, i.e. U = U{t). In SHE systems, the t-dependence of the 
unitary transformations naturally arises due to the acceleration of carriers in the presence of 
an electric field: the electron wave-vector (k) changes linearly in t, and consequently (B(k)) 
acquires a time- dependence. To incorporate the explicit time-dependence of the system 
quantum mechanically, we switch to the interaction picture.— In this picture, the original 
Hamiltonian ( TTTT) is split into two parts, H = Ho + Hi, where 

H = eE-f (20) 



governs the time evolution of the operators, and 



governs the time evolution of the states. In the usual sense, an operator A in the 
Schrodinger picture is transformed to the interaction picture (subscript I) as Ai(t) = 
e iHot/hj^ e -tHot/h an( j carr j es an explicit time-dependence by satisfying the Heisenberg re- 
lation, Ai = (z/i) _1 [v4/, Ho]. In particular, the momentum operator in the new picture is 
found to be pi(t) = p — eEt, i.e. with the expected linear time-dependence due to the elec- 
tric field. The state vectors \ip(t)) in the Schrodinger picture correspondingly transform as 
IV'i'(O) = e lHot ^ h \ip{t)) , and evolve according to the new "Schrodinger equation", 

n I (t)\Mt)) = m\^i(t)}, (22) 

where Hi(t) = e i ^ t / h n ie - iH o t / h . For the case of linear (e.g. Rashba) spin-orbit coupling, 
Hi(t) is evaluated to be 

^)^---(^)-^)^---^ <*> 

Higher order spin-orbit terms (~ k n , n > 2) generally lead to correspondingly higher order 
partial derivatives of the spin-orbit field in 7ii(t). The Hamiltonian (123!) governing the state 
vector evolution in the interaction picture is that of an electron subject to an explicitly time- 
dependent magnetic field, which we denote as Bit). Analogous to our previous treatment, we 
proceed to diagonalize the Schrodinger equation (1221) at time t, by applying a unitary rotation 
U(t) (defined as in Eq. ffl2l) but with 9 and in m carrying an explicit time- dependence). 
This transformation aligns the 5-axis to be parallel to the instantaneous magnetic field B(t), 
i.e. 

UtyH^U^t) = U(t)(iHd t )tf(t), 

" 2 

EL- 7<Tz \B{t)\ = ihB t + ihU{t)d t U\t), 
2m 

= e-Mo(t). (24) 

where e = ifidt is the energy operator. On the left-hand-side, the local transformation 
diagonalizes the time-dependent Zeeman term as required. On the right-hand-side, we obtain 
from the time-dependence of U a gauge field Ao(t) = —iU(t)d t W(t) related to the temporal 
evolution of the system. From the relations in Eqs. (fl4j) and (fl5]) . we can express the gauge 
field as Ao(t) = A t ■ <?, where A t = mxm = ^nxn + (A t ■ n)n. Thus, the term hAo(t) 
represents an additional Zeeman-like term, indicating the presence of an effective magnetic 
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field in the rotating frame. We elucidate the origin of this field in more detail below. The 
unitary transformation U(t) we invoked defines an instantaneous angular velocity u l = l (t) 
of the coordinates (in the laboratory frame, /) as it follows the time-dependent magnetic 
field. In the rotating frame r, this vector is given by uj t , where a ■ u r = U (a ■ oo l ) W. Since 
U = ^iUa ■ uj 1 [32], the Zeeman-like term hAo(t) in (j24p thus yields —fijlia ■ uj t ), which 
corresponds to an effective magnetic field — u r (omitting a scaling factor) in the rotating 
frame. This translates into an effective magnetic field B t = —uj 1 in the laboratory frame. If 
we now denote by n = n(t) the unit vector pointing along the direction of the magnetic field 
at time t, we have the equation of motion n = uj 1 x n. Performing a post cross product on 
both sides by n, one arrives at the expression for the angular velocity uj 1 = nxn + [uj 1 ■ n) n, 
or, in terms of the effective magnetic field, 

B t = kx n + (^Bt-njn. (25) 

Thus, the effective magnetic field arising from the gauge field Ao(t) of the unitary transfor- 
mation has a component along nxn and along n. Note that it does not have any component 
along n. As we noted previously, the unitary rotation matrix used by us (fP2l) is not unique. 
Specifically, different rotation matrices Ui, each specifying distinct angular velocities u\, 
can be used to align the reference 5-axis with the instantaneous magnetic field B(t); the 
freedom of choice here lies in determining the trajectory of the remaining x, y-axes. The 
second term on the right-hand-side of Eq. fl25l) reflects the particular choice of the gauge 
transformation U^. It is not an invariant of the gauge transformation (its magnitude being 
dependent on the particular gauge choice), and does not represent a physical field. How- 
ever, the first component nxn oi the effective magnetic field is invariant with respect to the 
gauge transformation, depending only on the time-dependence of the magnetic field B(t). 
This term can be understood to be a direct consequence of the time-dependent rotation of 
the axes.— ^ The same expression can be derived classically by directly comparing the spin 
vector in adjacent time frames^ — as a complement to the quantum derivation, the classical 
treatment is shown in detail in the Appendix. The nxn component represents a physical 
magnetic field which couples to the electron spins,— and, as we show below, is precisely 
the component which leads to the SHE in Rashba 2DEG systems. 
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III. ANALYSIS 



A. The intrinsic SHE due to Rashba SOC 

The Hamiltonian of conduction electrons in the Rashba system is given by Eq. fflHl) . 
Following our analysis above, the time- dependence of the effective Rashba field Br due to 
the electrons' motion in momentum space necessarily gives rise to a secondary component 
B± = n x n, where n = p~ x {py, — p x , 0) is the unit vector in the direction of Br. We 
assume a longitudinally applied electric field along the x-direction E = E x x, so that n = 
p _1 (0, eE x , 0). Because Br is strictly in-plane (i.e. it lies in the x,^-plane of the 2DEG), 
the term n x n represents an out-of-plane magnetic field component which is along the in- 
direction by convention. Next, we apply the adiabatic condition for the electron spins. In 
the ideal adiabatic limit, the magnetic field \Br\ is infinitely strong, so the spins always 
remain aligned to it as it varies with time. In reality, \Br\ is finite and there is a non-zero 
secondary component B±, and the relevant condition is \Br\ ^> \B±\, i.e. the electron spin 
is primarily aligned to Br, but with a small deviation along B±. In terms of the parameters 
of the Rashba system, the adiabatic condition reads as 

ak 2 \e\ ,_, 

> E x . (26) 

e 

Inserting typical values for the Rashba parameter a = 10~ n eVm and the Fermi wave- 
vector k = 10 s m -1 , we arrive at the condition E x <C~ 10 5 Vm™ 1 , which usually holds true 
in experiments. Assuming that the spin of electrons follow the direction of the net effective 
magnetic field, which is the sum of the spin-orbit field Br and the secondary component, 
the classical spin vector is given by 

where ± represents spin aligned parallel (+) or anti-parallel (— ) to the net field. To first 
order, the component of the spin along the i-direction is 

s z = ±-^^(fixn) -z, (28) 



where, to be consistent in units, the magnetic field in the denominator is defined in terms 
of its equivalent angular velocity. Note that in the convention above the + corresponds to 



12 



the ground state e_, whilst — corresponds to the eigenstate e+. In the adiabatic limit, the 



Eq. (1291) above describes a transverse separation of spins in the Rashba system. Let us 
consider the case for the ground state. Since the spin s z oc — p y , we find that electrons moving 
in the +i/-direction are polarized out-of-plane along the —z direction, whereas those moving 
in the —^-direction are polarized along +z. For the other eigenstate, the direction of the 
polarization is reversed, and hence there is a certain degree of canceling of the polarization 
if both eigenstates are present. However, at the Fermi level, there are more electrons in 
the ground state e_, giving rise to a net transverse spin separation and hence the SHE 
described in Ref. Q (see Eqs. (5) — (7) there). Summing over the Fermi surfaces of the 
two eigenstates yields an intrinsic spin-Hall (sH) conductivity of <t s h = jy/E x = — e/87r, 
where j* = h/4{s z ,v y } is the transverse spin-current. From our analysis above, we have 
clarified that the SHE in Rashba systems occurs as a result of an adiabatic process, in which 
electrons' spins become aligned to momentum-dependent magnetic fields that arise from the 
time-dependence of the system. The effect is therefore not due to the precessional behavior 
of spins, as is often stated in the literature. 

1. Berry's phase 

We alluded earlier to previous work which related the intrinsic spin-Hall conductivity in 
Rashba systems to the /c-space Berry phase of electrons through the Kubo formula.- It was 
found there that <t s h = eip±/8n 2 , where ip± is the Berry phase of electrons, 



magnitude of By, approaches that of Br, and applying this limit to Eq. (1281) . we obtain for 
the out-of-plane spin polarization 




eh 3 p y E x 
4ap 3 



(29) 



f 




(k) ■ dk 



2 



s 



(1 - cos6»)V k (/> • dk. 



(30) 



13 



The natural parameterization for the vector k is the time variable t, and rewriting the line 
integral above in terms of t we obtain 



Thus the Berry phase and hence the intrinsic spin-Hall conductivity of the Rashba system 
can be written equivalently in terms of the time component of the adiabatic gauge field, 



2. Effects of disorder 

Previous studies have shown that the intrinsic SHE in infinite Rashba systems vanishes 
in the presence of disorder.-^ Specifically, the vertex correction was shown to exactly cancel 
the intrinsic conductivity of e/87r even in the weak scattering limit. We provide a heuristic 
argument based on our analysis for the vanishing SHE. In the presence of disorder, the 
scattering provides a braking effect which cancels the acceleration of carriers on average in 
the steady state.— This implies that in the steady state we have (k) = 0, i.e. there is no 
net change in the momentum and thus the magnetic field component B± = n x n averages 
out to zero. Note, however, that this picture is an oversimplification,^ and that the SHE 
in Rashba systems does not vanish in general. For example, the SHE persists in finite-sized 
systems^^ 1 ^ and in the presence of spin-dependent impurities. 13 

B. The intrinsic SHE due to linear Dresselhaus SOC 

The case for the linear Dresselhaus spin-orbit coupling is also easily verified by our anal- 
ysis. The Dresselhaus spin-orbit Hamiltonian is given by 




(31) 



it). 




(32) 



where (3 is the Dresselhaus SOC strength and Br> is the effective Dresselhaus SOC field. 
Here we have n = p~ 1 (p x , —p y ,0), and we find that {n x n) z = +eE x p y /p 2 . Consequently, 
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the out-of-plane spin polarization s z has the same magnitude but opposite sign compared to 
the Rashba SOC case. This is in agreement with previous theoretical studies^ which predicts 
the spin-Hall conductivity in this system to be <t s h = e/871. 

IV. DISCUSSIONS 

Having established the Rashba SHE as an adiabatic effect, we now identify two common 
traits of the two intrinsic SHEs; adiabaticity and time-dependence. For the Rashba system, 
we found that the time- dependent spin-orbit field B R {t) = \B R {t)\n is always accompanied 
by an additional effective magnetic field, B± = n x n. This correction to the magnetic field 
results in a net field By, = Br + B± which is different to Br. Considering the adiabatic 
limit, where \Br\ ^> we recovered exactly the results of Sinova et al. describing the 

SHE in the Rashba system. The field B± was shown to be described quantum mechanically 
by a gauge field in time space. In the Luttinger system, the adiabatic assumption results 
in a non-trivial momentum space curvature which enters the equation of motion as the 
spin-dependent anomalous velocity component ({TO]) . 

Thus, at first glance it appears that the two effects are rather independent phenomena. 
However, an interesting duality exists between the two effects. In the former Luttinger case, 
a spin-dependent anomalous velocity pushes opposite spin species to opposite lateral sides of 
the sample. The magnetic field responsible for this effect is the Berry curvature defined by 
the /c-space gauge field. On the other hand, in the Rashba system, a momentum- dependent 
magnetic field polarizes electrons along opposite directions out-of-the-plane depending on 
their transverse propagation direction. The magnetic field responsible for this effect is 
defined by the i-space gauge field. Given this duality, it would be tempting to ask whether 
there is any underlying relation between the two pictures. 

We proceed to consolidate the link between the two effects by investigating the connection 
between the anomalous velocity due to the Berry curvature, and the presence of the B± term. 
In this endeavour, we employ the reciprocal space analogue of the analysis by Aharanov and 
Stern^ of the origin of the Berry's curvature in real space. We consider again the general 
spin-orbit Hamiltonian in Eq. ( fTTT) . The velocity along the i-th coordinate is given by 
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Hamilton's equation: 

in 

Pi dB(k) _ 

= 7^ o- (33) 

m opi 

When the magnetic field is time-dependent, the spins see an additional magnetic field B±. 
Assuming that spins align to B% — B + B±, the spin vector, to first order, is given by 
a = By, /I By, |. Writing the spin-orbit field B = \B\n, the partial derivative in Eq. (133]) can 
be expanded into its magnitude and directional parts as (d\B\/dpi)n + \B\dn/dpi. Taking 
the adiabatic limit |£?£|/|J3| — > 1, the second term in the velocity expression becomes 

d\B\ h / \ dn 
Vi = ^ - 2 (" X ") ■ Wi (M) 

Writing n = kjdn/dkj, where the summation over j is implicit, and rearranging the terms 
we then get 

d\B\ kj f dn dn\ _ 



dpi 2 \dki dkj y 
The first term in the above equation represents a velocity term that is due to the inhomo- 
geneity of the spin-orbit field B in momentum space, i.e. it is the reciprocal space analogue 
of the Stern-Gerlach force. Remarkably, the second term in Eq. (|35l) is the anomalous ve- 
locity of electrons due to Berry's curvature in /c-space. This becomes clearer when written 
in terms of the magnetic field vector B = \B\n, 

v—- - k g (9B dB 

= -eijkkjVtkik). (36) 

We find that this is exactly the anomalous velocity component in Eq. (Q, which arises 
from the Berry's curvature in fc-space. By letting the spin vector to be aligned anti-parallel 
to the net magnetic field, a = —By/\By\ in Eq. fl33l) . the correct spin-dependence of the 
anomalous velocity can be obtained. Thus, we have shown that the anomalous velocity due 
to the Berry curvature actually arises because of the B± magnetic field component, which 
in turn is related to the time component of the gauge field. This gauge field component 
therefore plays an equally important role in the SHE in the Luttinger system, as it does in 
the Rashba one, and acts as the unifying bridge between the two effects. 
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V. SUMMARY 



The primary motivation for this paper is to establish the link between the two intrinsic 



SHEs reported in Refs. 



au|) which has not been clarified hitherto. We first considered the 
intrinsic SHE in the Luttinger system, which is driven by the spin-dependent anomalous 
velocity due to the non-trivial curvature of momentum space. However, this theoretical 
picture is not applicable in the planar Rashba system. Instead, the SHE in the Rashba 
system was shown to arise from spins acquiring a component (in the adiabatic sense) along 
an additional effective magnetic field B±, arising from the time-dependence of the system. 
This field component was shown to be described by a gauge field in time space. Finally, we 
showed that in the adiabatic limit, B± is also the origin of the anomalous velocity due to the 
momentum space Berry curvature. Thus, we conclude that the intrinsic SHEs in the two 
systems are simply different manifestations of B±, and that this term provides a unifying 
link between the two effects. 



APPENDIX: CLASSICAL DERIVATION OF EFFECTIVE MAGNETIC FIELD 
COMPONENT, k x n 

Consider the dynamics of the spin vector s(t) in a time-dependent magnetic field B(t), 

m=g(s(t)xB(t)), (A.l) 

where g is the coupling factor. To solve the above equation, we freeze the time-dependence 
by transforming to a rotated coordinate frame at each point in time, such that the 5-axis is 
aligned with the magnetic field. A spin vector s defined relative to the coordinate frame at 
time t, is expressed as the vector s ' = s + sx uj(t)dt in the coordinate frame at time t + dt, 
where u(t) is the generator of infinitesimal rotations (see Fig. [1]). The choice of uj(t) is not 
unique; however, specifically choosing i2{t) = zx z where z is the unit vector n = B/\B\ as 
seen in the rotated frame, coincides with the parallel transport of the coordinate frames.— ^ 
Suppose we have a vector representing the spin, sit), in the rotated frame at time t. At 
time t + dt, this vector becomes [relative to frame t + dt] s(t + dt) + s(t + dt) x uo{t)dt where 
s(t + dt) « s(t) + g{s(t) x \B\z)dt. For infinitesimally small dt, we may write s(t + dt) [in 
frame t] m s(t + dt) + s(t + dt) x uj(t)dt [in frame t + dt]. The right-hand-side of the resulting 
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FIG. 1: (Color online) (left) The classical spin vector s(t) precesses about the magnetic field which 
is along the z direction at some instant t. Because of the time-dependence of the magnetic field, 
the spin is also subject to a rotation about u)(t) = z x z (see text) which transforms it from the 
frame at time t (left) to the frame at time t + dt (right). The cD(t) acts as an additional magnetic 
field which governs the overall spin dynamics. 

equation becomes s(t) + g(sx \B\z)dt + s x uo{t)dt + 0(dt 2 ). Rearranging, and taking the 
limit dt — > 0, we have 

Hm s(t + dt)-s(t) = g = ( gx \B(t)\z) +sx u{t), 
dt^o dt \ / w 

= sx (g\B(t)\z '+ 'z x zj . (A.2) 

Therefore, as seen in the laboratory frame, there is an additional, effective magnetic field 
_B e ff. = B{t) + g~ x n x n. 
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36 This argument was stated previously in the context of spin precession in 14]: in the presence 
of impurities, the scattering scrambles the spin precession sufficiently such that no net SHE 
results. 
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